Examination of the Potential of Terrestrial Laser Scanning and Structure-from-Motion Photogrammetry for Rapid Nondestructive Field Measurement of Grass Biomass

Above ground biomass (AGB) is a parameter commonly used for assessment of grassland systems. Destructive AGB measurements, although accurate, are time consuming and are not easily undertaken on a repeat basis or over large areas. Structure-from-Motion (SfM) photogrammetry and Terrestrial Laser Scanning (TLS) are two technologies that have the potential to yield precise 3D structural measurements of vegetation quite rapidly. Recent advances have led to the successful application of TLS and SfM in woody biomass estimation, but application in natural grassland systems remains largely untested. The potential of these techniques for AGB estimation is examined considering 11 grass plots with a range of biomass in South Dakota, USA. Volume metrics extracted from the TLS and SfM 3D point clouds, and also conventional disc pasture meter settling heights, were compared to destructively harvested AGB total (grass and litter) and AGB grass plot measurements. Although the disc pasture meter was the most rapid method, it was less effective in AGB estimation (AGBgrass r2 = 0.42, AGBtotal r2 = 0.32) than the TLS (AGBgrass r2 = 0.46, AGBtotal r2 = 0.57) or SfM (AGBgrass r2 = 0.54, AGBtotal r2 = 0.72) which both demonstrated their utility for rapid AGB estimation of grass systems

Comparisons of Temporal and Spatial Trends in the Spatially Complete Global Spectral Surface Albedos Products

Five years of spatially complete snow-free land surface albedo data have been prepared using high quality white-sky and black-sky land surface albedo observations (MOD43B3) from the MODIS instrument aboard NASA's Terra satellite platform. The data were generated using an updated ecosystem-dependent temporal interpolation technique. In this paper, we describe the refinements in the technique and the creation of a spatially complete snow-free five-year aggregate climatology product. This paper also describes an error analysis of the interpolation technique. The filled albedo products are examined through comparisons of temporal and spatial trends for pixels that have been filled versus pixels that have been retained in the original MOD43B3 values. The variability in the trends showcase how the filling technique maintains the pixel-level spatial, spectral, and temporal integrity of the MOD43B3 data. These comparisons are made for both a single-year of filled data, year 2002, and for the five-year aggregate climatology product

Assessment of Cloud Screening with Apparent Surface Reflectance in Support of the ICESat-2 Mission

The separation of cloud and clear scenes is usually one of the first steps in satellite data analysis. Before deriving a geophysical product, almost every satellite mission requires a cloud mask to label a scene as either clear or cloudy through a cloud detection procedure. For clear scenes, products such as surface properties may be retrieved; for cloudy scenes, scientist can focus on studying the cloud properties. Hence the quality of cloud detection directly affects the quality of most satellite operational and research products. This is certainly true for the Ice, Cloud, and land Elevation Satellite-2 (lCESat-2), which is the successor to the ICESat-l. As a top priority mission, ICESat-2 will continue to provide measurements of ice sheets and sea ice elevation on a global scale. Studies have shown that clouds can significantly affect the accuracy of the retrieved results. For example, some of the photons (a photon is a basic unit of light) in the laser beam will be scattered by cloud particles on its way. So instead of traveling in a straight line, these photons are scattered sideways and have traveled a longer path. This will result in biases in ice sheet elevation measurements. Hence cloud screening must be done and be done accurately before the retrievals

Quantifying the effects of freeze-thaw transitions and snowpack melt on land surface albedo and energy exchange over Alaska and Western Canada

Variations in land surface albedo and snow-cover strongly impact the global biosphere, particularly through the snow-albedo feedback on climate. The seasonal freeze-thaw (FT) transition is coupled with snowpack melt dynamics and strongly impacts surface water mobility and the energy budget in the northern (≥45°N) arctic and boreal region (ABR). However, understanding of the regional variation in snowmelt and its effect on the surface energy budget are limited due to sparse in situ measurements of these processes and environmental constraints on effective monitoring within the ABR. In this study, we combined synergistic observations from overlapping satellite optical-infrared and microwave sensor records to quantify the regional patterns and seasonal progression in wet snow conditions during the spring snowmelt and autumn snow accumulation periods across Alaska and western Canada. The integrated satellite record included daily landscape FT status from AMSR microwave brightness temperature retrievals; and snow-cover extent, black sky albedo and net shortwave solar radiation (R snet) derived from MODIS and AVHRR observations. The integrated satellite records were analyzed with in situ surface air temperature and humidity observations from regional weather stations over a two-year study period (2015–2016) overlapping with the NASA ABoVE (Arctic Boreal Vulnerability Experiment). Our results show a large (79%) mean decline in land surface albedo between dry snow and snow-free conditions during the spring (March–June) and autumn (August–November) transition periods. Onset of diurnal thawing and refreezing of the surface snow layer and associated wet snow conditions in spring contributed to an approximate 25% decrease in snow cover albedo that extended over a seven to 21 week snowpack depletion period. The lower wet snow albedo enhances R snet by approximately 74% (9–10 MJ m−2 d−1) relative to dry snow conditions, reinforcing snowmelt and surface warming, and contributing to growing season onset and activation of biological and hydrological processes in the ABR. These results contribute to better understanding of snow albedo feedbacks to Arctic amplification, and the representation of these processes in global Earth system models

An Approach for the Long-Term 30-m Land Surface Snow-Free Albedo Retrieval from Historic Landsat Surface Reflectance and MODIS-based A Priori Anisotropy Knowledge

Land surface albedo has been recognized by the Global Terrestrial Observing System (GTOS) as an essential climate variable crucial for accurate modeling and monitoring of the Earth's radiative budget. While global climate studies can leverage albedo datasets from MODIS, VIIRS, and other coarse-resolution sensors, many applications in heterogeneous environments can benefit from higher-resolution albedo products derived from Landsat. We previously developed a "MODIS-concurrent" approach for the 30-meter albedo estimation which relied on combining post-2000 Landsat data with MODIS Bidirectional Reflectance Distribution Function (BRDF) information. Here we present a "pre-MODIS era" approach to extend 30-m surface albedo generation in time back to the 1980s, through an a priori anisotropy Look-Up Table (LUT) built up from the high quality MCD43A BRDF estimates over representative homogenous regions. Each entry in the LUT reflects a unique combination of land cover, seasonality, terrain information, disturbance age and type, and Landsat optical spectral bands. An initial conceptual LUT was created for the Pacific Northwest (PNW) of the United States and provides BRDF shapes estimated from MODIS observations for undisturbed and disturbed surface types (including recovery trajectories of burned areas and non-fire disturbances). By accepting the assumption of a generally invariant BRDF shape for similar land surface structures as a priori information, spectral white-sky and black-sky albedos are derived through albedo-to-nadir reflectance ratios as a bridge between the Landsat and MODIS scale. A further narrow-to-broadband conversion based on radiative transfer simulations is adopted to produce broadband albedos at visible, near infrared, and shortwave regimes.We evaluate the accuracy of resultant Landsat albedo using available field measurements at forested AmeriFlux stations in the PNW region, and examine the consistency of the surface albedo generated by this approach respectively with that from the "concurrent" approach and the coincident MODIS operational surface albedo products. Using the tower measurements as reference, the derived Landsat 30-m snow-free shortwave broadband albedo yields an absolute accuracy of 0.02 with a root mean square error less than 0.016 and a bias of no more than 0.007. A further cross-comparison over individual scenes shows that the retrieved white sky shortwave albedo from the "pre-MODIS era" LUT approach is highly consistent (R(exp 2) = 0.988, the scene-averaged low RMSE = 0.009 and bias = 0.005) with that generated by the earlier "concurrent" approach. The Landsat albedo also exhibits more detailed landscape texture and a wider dynamic range of albedo values than the coincident 500-m MODIS operational products (MCD43A3), especially in the heterogeneous regions. Collectively, the "pre-MODIS" LUT and "concurrent" approaches provide a practical way to retrieve long-term Landsat albedo from the historic Landsat archives as far back as the 1980s, as well as the current Landsat-8 mission, and thus support investigations into the evolution of the albedo of terrestrial biomes at fine resolution

Evaluation of Surface and Near-Surface Melt Characteristics on the Greenland Ice Sheet using MODIS and QuikSCAT Data

The Greenland Ice Sheet has been the focus of much attention recently because of increasing melt in response to regional climate warming. To improve our ability to measure surface melt, we use remote-sensing data products to study surface and near-surface melt characteristics of the Greenland Ice Sheet for the 2007 melt season when record melt extent and runoff occurred. Moderate Resolution Imaging Spectroradiometer (MODIS) daily land-surface temperature (LST), MODIS daily snow albedo, and a special diurnal melt product derived from QuikSCAT (QS) scatterometer data, are all effective in measuring the evolution of melt on the ice sheet. These daily products, produced from different parts of the electromagnetic spectrum, are sensitive to different geophysical features, though QS- and MODIS-derived melt generally show excellent correspondence when surface melt is present on the ice sheet. Values derived from the daily MODIS snow albedo product drop in response to melt, and change with apparent grain-size changes. For the 2007 melt season, the QS and MODIS LST products detect 862,769 square kilometers and 766,184 square kilometers of melt, respectively. The QS product detects about 11% greater melt extent than is detected by the MODIS LST product probably because QS is more sensitive to surface melt, and can detect subsurface melt. The consistency of the response of the different products demonstrates unequivocally that physically-meaningful melt/freeze boundaries can be detected. We have demonstrated that these products, used together, can improve the precision in mapping surface and near-surface melt extent on the Greenland Ice Sheet

Mid-summer snow-free albedo across the Arctic tundra was mostly stable or increased over the past two decades

Arctic vegetation changes, such as increasing shrub-cover, are expected to accelerate climate warming through increased absorption of incoming radiation and corresponding decrease in summer shortwave albedo. Here we analyze mid-summer shortwave land-surface albedo and its change across the pan-Arctic region based on MODIS satellite observations over the past two decades (2000-2021). In contrast to expectations, we show that terrestrial mid-summer shortwave albedo has not significantly changed in 82% of the pan-Arctic region, while 14% show an increase and 4% a decrease. The total median significant change was 0.014 cumulative over the past 22 years. By analyzing the visible and near-/shortwave-infrared range separately, we demonstrate that the slight increase arises from an albedo increase in the near-/shortwave infrared domain while being partly compensated by a decrease in visible albedo. A similar response was found across different tundra vegetation types. We argue that this increase in reflectance is typical with increasing biomass as a result of increased multiple reflection in the canopy. However, CMIP6 global climate model albedo predictions showed the opposite sign and different spatial patterns of snow-free summer albedo change compared to satellite-derived results. We suggest that a more sophisticated vegetation parametrization might reduce this discrepancy, and provide albedo estimates per vegetation type

Variability in Surface BRDF at Different Spatial Scales (30 m-500 m) Over a Mixed Agricultural Landscape as Retrieved from Airborne and Satellite Spectral Measurements

Over the past decade, the role of multiangle remote sensing has been central to the development of algorithms for the retrieval of global land surface properties including models of the bidirectional reflectance distribution function (BRDF), albedo, land cover/dynamics, burned area extent, as well as other key surface biophysical quantities represented by the anisotropic reflectance characteristics of vegetation. In this study, a new retrieval strategy for fine-to-moderate resolution multiangle observations was developed, based on the operational sequence used to retrieve the Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 5 reflectance and BRDF/albedo products. The algorithm makes use of a semiempirical kernel-driven bidirectional reflectance model to provide estimates of intrinsic albedo (i.e., directional-hemispherical reflectance and bihemispherical reflectance), model parameters describing the BRDF, and extensive quality assurance information. The new retrieval strategy was applied to NASA's Cloud Absorption Radiometer (CAR) data acquired during the 2007 Cloud and Land Surface Interaction Campaign (CLASIC) over the well-instrumented Atmospheric Radiation Measurement Program (ARM) Southern Great Plains (SGP) Cloud and Radiation Testbed (CART) site in Oklahoma, USA. For the case analyzed, we obtained approx.1.6 million individual surface bidirectional reflectance factor (BRF) retrievals, from nadir to 75 off-nadir, and at spatial resolutions ranging from 3 m - 500 m. This unique dataset was used to examine the interaction of the spatial and angular characteristics of a mixed agricultural landscape; and provided the basis for detailed assessments of: (1) the use of a priori knowledge in kernel-driven BRDF model inversions; (2) the interaction between surface reflectance anisotropy and instrument spatial resolution; and (3) the uncertain ties that arise when sub-pixel differences in the BRDF are aggregated to a moderate resolution satellite pixel. Results offer empirical evidence concerning the influence of scale and spatial heterogeneity in kernel-driven BRDF models; providing potential new insights into the behavior and characteristics of different surface radiative properties related to land/use cover change and vegetation structure

Seasonal variation of source contributions to eddy-covariance CO 2 measurements in a mixed hardwood-conifer forest

Net ecosystem exchange (NEE) measurements using the eddy covariance technique have been widely used for calibration and evaluation of carbon flux estimates from terrestrial ecosystem models as well as for remote sensing-based estimates across various spatial and temporal scales. Therefore, it is vital to fully understand the land surface characteristics within the area contributing to these flux measurements (i.e. source area) when upscaling plot-scale tower measurements to regional-scale ecosystem estimates, especially in heterogeneous landscapes, such as mixed forests. We estimated the source area of a flux tower at a mixed forest (Harvard Forest in US) using a footprint model, and analyzed the spatial representativeness of the source area for the vegetation characteristics (density variation and magnitude) within the surrounding 1- and 1.5-km grid cells during two decades (1993–2011). Semi-variogram and window size analyses using 19 years of Landsat-retrieved enhanced vegetation index (EVI) confirmed that spatial heterogeneity within the 1-km grid cell has been gradually increasing for leaf-on periods. The overall prevailing source areas lay toward the southwest, yet there were considerable variations in the extents and the directions of the source areas. The source areas generally cover a large enough area to adequately represent the vegetation density magnitude and variation during both daytime and nighttime. We show that the variation in the daytime NEE during peak growing season should be more attributed to variations in the deciduous forest contribution within the source areas rather than the vegetation density. This study highlights the importance of taking account of the land cover variation within the source areas into gap-filling and upscaling procedures